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ABSTRACT 

We present bright galaxy number counts measured with linear detectors in the 
B, V, I, and K bands in two fields covering nearly 10 square degrees. All of our 
measurements are consistent with passive evolution models, and do not confirm the 
steep slope measured in other surveys at bright magnitudes. Throughout the range 
16 < B < 19, our B— band counts are consistent with the "high normalization" 
models proposed to reduce the faint blue galaxy problem. Our A'— band counts agree 
with previous measurements, and have reached a fair sample of the universe in the 
magnitude range where evolution and K-corrections are well understood. 
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1 INTRODUCTION 

Observational studies of galaxy formation and evolution 
have advanced at an unprecedented pace in recent years. 
Two developments have played a key role: CCD imagery 
to very faint limits and the ability to measure redshifts 
for large samples at increasingly faint magnitudes. Progress 
to date has relied primarily on optical data, but it has 
been clear for some time that near-infrared observations 
are fundamental. Samples selected according to A— band 
flux are superior to the traditional B— selected samples in 
at least three respects: (i) In the infrared, K-corrections 
due to the redshift of the spectral energy distribution are 
smooth, well-understood, and nearly independent of Hub- 
ble type; the expected luminosity evolution is also smooth. 
(ii) At high redshift, the observer's near-infrared samples 
the well-understood rest frame optical (dominated by long- 
lived, near-solar mass stars), while the optical band samples 
the poorly-understood rest frame ultraviolet (dominated by 
short-lived massive stars.) (Hi) Since near-solar mass stars 
make up the bulk of a galaxy, the absolute A magnitude is 
a measure of the visible mass in a galaxy. 

Deep photometric surveys of small areas measuring the 
number counts and colours of field galaxies have reached as 
faint as A = 24, (Gardner, Cowie & Wainscoat 1993; Cowie 
et al. 1994; Djorgovski et al. 1995,) and show significant 
amounts of galaxy evolution at high redshift. However, the 



interpretation of faint galaxy data (either counts or redshift 
distributions) hinges on an accurate statistical description 
of the local population of galaxies. Photometric and spec- 
troscopic surveys of bright galaxies, covering an area large 
enough to average over the effects of large-scale structure, 
are required in order to obtain bright galaxy counts and 
colour distributions, measurements of the galaxy-galaxy cor- 
relation function, and the local luminosity function. In the 
A band, the small size of infrared detectors and the corre- 
sponding small field of view available, has made it difficult in 
the past to image large areas. In the optical, much work has 
been done using photographic plates, but their non-linearity 
can introduce possible systematics in the photometry (Met- 
calfe, Fong & Shanks 1995.) 

We have imaged nearly 10 square degrees in two fields 
at high galactic latitude using a NICMOS3 detector in the 
near-infrared A band, and a CCD camera in the optical 
B, V, and / bands. Here, we present the galaxy number 
counts. In a companion paper (Baugh et al. 1996, Paper 
II) we present the galaxy correlation function. Our A < 
15 photometric catalog has been used to select galaxies for 
spectroscopic follow-up, and future papers in this series will 
present the A— band galaxy luminosity function, the galaxy 
redshift and colour distributions, and a discussion of the star 
counts and colour distribution. One of our fields, centered 
on the north ecliptic pole, is ideally situated for viewing by 
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satellites in polar orbits, and has been the subject of deep 
IRAS and ROSAT observations. 



2 THE DATA 

We have imaged 9.84 square degrees in the A"— band with 
a 256 2 HgCdTe NICM0S3 detector, and in the B, V, and 
I bands with a 2048 2 CCD camera. The if— band obser- 
vations were made in 1994 June with the IRIM camera on 
the Kitt Peak National Observatory 1.3m telescope. On this 
telescope, the IRIM camera has 1.96" pixels and an 8.36' 
field of view. Each point was observed with at least two 
60sec exposures, reaching a 5a galaxy detection depth of 
K = 15.6 in a 10" circular aperture. The 3a surface bright- 
ness limit of our images is K = 18.5 mag arcsec~ 2 . The 
B-, V— , and /—band observations were made in 1995 June 
with the T2KA camera on the KPNO 0.9m telescope. On 
this telescope, the T2KA CCD has 0.68" pixels with a 23.2' 
field of view. Each point in each filter was observed with a 
300sec exposure, and the images reach a 5a detection depth 
of B = 21.1, V = 20.9, and I = 19.6 in a 10" circular 
aperture. The 3a surface brightness limit is 24.0, 23.8, and 
22.5 mag arcsec~ for B, V and I respectively. The location 
of the two fields were selected randomly, that is, without re- 
gard to the presence or absence of any known objects. The 
field centers are at RA 14 h 15 m , Dec +00 and RA 18 h 00 m , 
Dec +66; galactic latitudes +55 and +30 respectively. One 
of our fields has a nearby rich galaxy cluster within it, and to 
avoid biasing the galaxy counts, we have removed all objects 
within 1 degree radius of the central galaxy of this cluster. 
Thus the effective area for the counts presented here is 8.54 
square degrees. 

Image reduction will be described in detail else- 
where (Gardner et al, in preparation), but mainly followed 
the techniques described by Gardner (1995). Briefly, the 
if— band images were dark-subtracted, flattened with dome- 
flats, median sky-subtracted, and large-scale gradients were 
removed with a large (2') median filter. The optical images 
were bias subtracted, and flattened with twilight flats. In 
addition, the i— band images were flattened with median 
sky flats. Object identification on the optical images was 
done with the SExtractor program (Bertin & Arnouts 1996), 
using a 3a threshold. Deblending was done with a multi- 
thresholding approach, and integral pixels were assigned to 
each object. Magnitudes are the SExtractor mag-best, which 
is the flux in an elliptical aperture at 2.5 times the Kron 
(1980) radius along the semi-major and semi-minor axes for 
isolated objects, and a corrected-isophotal magnitude when 
deblending was required, if— band photometry was done in 
10" apertures on all identified I— band objects. The aperture 
magnitudes were deblended by assigning pixels by hand, and 
were corrected to total magnitudes using the curve of growth 
measured in the I band. Colours for all objects were mea- 
sured within 5" circular apertures in the optical, and within 
6" square apertures in I — K . 

Star/galaxy separation is critical at these magnitudes, 
and was carried out using a combination of morphology on 
the optical images, and the colour criterion discussed by 
Gardner (1995). The seeing of the optical images varied 
through the observing run between two and three pixels, or 
1.3" < fwhm < 2.0". All identified galaxies with K < 15, 



Table 1. The galaxy number counts 



Mag 


Raw N 


log(N) 


&high 




10.25 


1 


-0.630 


0.519 


0.762 


10.75 


1 


-0.630 


0.519 


0.762 


11.25 


4 


-0.028 


0.253 


0.283 


11.75 


13 


0.484 


0.134 


0.139 


12.25 


22 


0.712 


0.101 


0.103 


12.75 


33 


0.888 


0.082 


0.083 


13.25 


66 


1.189 


0.057 


0.057 


13.75 


138 


1.510 


0.039 


0.039 


14.25 


273 


1.806 


0.027 


0.027 


14.75 


642 


2.177 


0.018 


0.018 


15.25 


1290 


2.480 


0.012 


0.012 


15.75 


2609 


2.786 


0.009 


0.009 


12.25 


2 


-0.329 


0.365 


0.451 


12.75 


1 


-0.630 


0.519 


0.762 


13.25 


6 


0.148 


0.203 


0.219 


13.75 


11 


0.411 


0.147 


0.153 


14.25 


23 


0.731 


0.099 


0.101 


14.75 


26 


0.785 


0.093 


0.094 


15.25 


63 


1.169 


0.058 


0.058 


15.75 


95 


1.347 


0.047 


0.047 


16.25 


198 


1.666 


0.032 


0.032 


16.75 


398 


1.970 


0.022 


0.022 


17.25 


644 


2.179 


0.017 


0.018 


17.75 


1190 


2.445 


0.013 


0.013 


12.25 


1 


-0.630 


0.519 


0.762 


13.25 


1 


-0.630 


0.519 


0.762 


13.75 


1 


-0.630 


0.519 


0.762 


14.25 


5 


0.069 


0.224 


0.246 


14.75 


6 


0.148 


0.203 


0.219 


15.25 


14 


0.516 


0.129 


0.133 


15.75 


25 


0.768 


0.095 


0.096 


16.25 


48 


1.051 


0.067 


0.067 


16.75 


83 


1.289 


0.050 


0.050 


17.25 


142 


1.522 


0.038 


0.038 


17.75 


285 


1.824 


0.026 


0.027 


18.25 


454 


2.027 


0.021 


0.021 


18.75 


766 


2.254 


0.016 


0.016 


12.25 


1 


-0.630 


0.519 


0.762 


14.25 


3 


-0.153 


0.295 


0.341 


14.75 


2 


-0.329 


0.365 


0.451 


15.25 


6 


0.148 


0.203 


0.219 


15.75 


6 


0.148 


0.203 


0.219 


16.25 


20 


0.671 


0.106 


0.109 


16.75 


26 


0.785 


0.093 


0.094 


17.25 


56 


1.118 


0.062 


0.062 


17.75 


117 


1.438 


0.042 


0.042 


18.25 


188 


1.644 


0.033 


0.033 


18.75 


323 


1.879 


0.025 


0.025 


19.25 


509 


2.076 


0.020 


0.020 


19.75 


901 


2.324 


0.015 


0.015 



The raw number of galaxies in the 8.54 square degree area, and 
the log(N/mag/deg 2 ). The high and low Poissonian errors are 
taken from the calculations of Gehrels (1986). 

i < 18, V < 19 or B < 20 were confirmed by eye. The 
median Kron (1980) n radius in the faintest magnitude bin 
on the optical images was typically 2.5" for stars and 3.4" 
for galaxies. The I — K colour, in combination with the 
B — I colour, is a good indicator of star/galaxy separation 
for all but the bluest objects. We found no large population 
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Figure 1. The K— band galaxy counts. The Euclidean slope, 
dlog(n) I 'dm = 0.6, has been subtracted in order to expand the 
ordinate. The models plotted are no-evolution (dotted line) and 
passive-evolution (solid line), in an open universe (the two middle 
lines at faint magnitudes) and in a closed universe. The dashed 
line is the no-evolution, local underdensity model of Huang et al. 
(1996). 



of compact objects with the colours of galaxies, nor did we 
find a large population of extended objects with the colours 
of stars. In the range 15 < if < 16, star/galaxy separation 
was carried out on the basis of colour alone, using V — I vs 
I — K for the objects not detected in B. 

The galaxy number counts are presented in Table |l[ 
and the if— band counts are plotted in Figure |l|. To ex- 
pand the ordinate, we have subtracted the Euclidean slope 
dlog(n) I dm — 0.6. Alongside our data, we plot other exist- 
ing bright if— band galaxy counts. For comparison we also 
show the predictions of a simple model, based upon the for- 
mulation of Yoshii & Takahara (1988), modified to include 
rest-frame and evolved spectral energy distributions from 
the GISSEL models (Bruzual & Chariot 1993; 1996 in prepa- 
ration.) The model is similar to that used in Gardner (1996), 
except that we have used the revised version of the GIS- 
SEL models. For the current purposes, the main difference 
between the two versions is that rest frame optical-near in- 
frared colours are redder. The solid lines include this passive 
evolution, while the dotted lines are no-evolution models, i.e. 
models that include only the cosmological geometry and K- 
corrections. The dashed line is a "local-underdensity" model 
proposed by Huang et al. (1996) and is discussed below. To 
constuct the models we adopted the bj type-independent 
luminosity function of Loveday et al. (1992), converted to 
type-dependent luminosity functions in other filters through 
rest frame colours. The normalization of the models was de- 
termined with a least-squares fit of the passive-evolution flat 
universe model to our data. 

Figure ^ shows our J— band galaxy counts, again with 
the Euclidean slope subtracted, together with other existing 



Figure 2. The /—band number counts. The models plotted are 
as in Figure All counts have been converted to as discussed 
in the text. The Medium Deep Survey counts of Casertano et al. 
(1995) show an excess over our counts and over those of the CFRS 
(Le Fevre et al. 1995). 



counts converted to the Kron-Cousins / band which we used. 
Figure ^ shows our V— band galaxy counts and Figure ^ our 
B— band galaxy counts. 



3 DISCUSSION 

The AT— band galaxy counts presented in Table |l| and Fig- 
ure |l| show remarkable consistency with the counts of Huang 
et al. (1996), which were also based upon if— band observa- 
tions of approximately 10 square degrees. While one of us 
(JPG) is also a co-author of that paper, the two data sets 
were collected, reduced and analyzed independently, and 
there is no overlap in area between the two surveys. The con- 
sistency of the two results indicates that both surveys rep- 
resent a fair sample of the universe, so the bright if— band 
galaxy counts and their normalization are no longer a sub- 
ject of debate. 

Our counts in the region 13 < if < 15 are within lcr 
of the Huang et al. (1996) counts in each magnitude bin. 
Nevertheless, we measure a shallower slope. The slope of 
our counts is 0.627 ±0.010, consistent with the flat-universe 
passive evolution model plotted in Figure [l]. Thus our data 
do not support the conclusion of Huang et al. (1996) that 
a model with a steep slope is required to fit the bright 
if —band number counts. The value of the normalization, 
4>* = 1.50 x 10" 2 /i 3 Mpc -3 , (for ii = WOh km s _1 Mpc' 1 ,) 
inferred from our counts is higher than the value in the 
Loveday et al. (1992) luminosity function converted to the 
if— band. This comparison, however, is uncertain because 
the transformation from the 6j-band luminosity function to 
a if— band luminosity function is very model dependent. 
The model if* depends strongly on the assumed rest-frame 
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colours and the 3 parameters of the Schechter (1976) lu- 
minosity function are strongly correlated. Neither of the 
two existing measurements of the K— band luminosity func- 
tion surveyed enough galaxies to accurately constrain <f>* 
(Mobasher, Sharpies & Ellis 1993; Glazebrook et al. 1995). 
In the absence of an accurate measurement of the joint 
[B, K] luminosity function, it is difficult to determine the 
consistency of the normalization of number count models 
in different filters. For this reason, we have normalized the 
models plotted in each of the figures to our data. 

To fit the bright K— band number counts without pas- 
sive evolution, Huang et al. (1996) have proposed that the 
local universe is underdense by a factor of 2. They have 
constructed a heuristic model in which this underdensity 
smoothly dampens out by z m 0.2. This model is unphysical 
and, as they stress, the number counts alone do not contain 
enough information to constrain its specifics. Maddox et al. 
(1990) measured a steep slope in the bright B— band num- 
ber counts and a local underdensity has also been proposed 
to explain this slope (Shanks 1990; Metcalfe et al. 1991). We 
plot the model of Huang et al. (1996) in Figure [l] (with their 
normalization) as a dashed line. This model does not fit our 
galaxy counts. 

The /—band galaxy counts are presented in Table [l] 
and Figure ^. We used the Kitt Peak standard I filter, 
corrected to the Kron-Cousins standard stars of Landolt, 
(1983; 1992). Also plotted in Figure |^ is a compilation of 
7— band counts, converted to Kron-Cousins Ike, taken from 
the Canada-France Redshift Survey (Le Fevre et al. 1995, 
hereafter CFRS) and the WF/PC Medium Deep Survey 
counts of Casertano et al. (1995), (hereafter MDS). The lat- 
ter were converted using I kc = /78s +0.16* (V555 — /78s) and 
(V 555 - Jres) = 1-6 (Edvardsson & Bell 1989). While there 
are no other measurements of the J— band number counts as 
bright as ours, the faint end of our counts are consistent with 
the CFRS counts. They are not, however, consistent with 
the MDS counts, which are a factor of 1.5 higher. The MDS 
measurements are based upon pre-refurbishment WF/PC 
images. Casertano et al. (1995) found that the median half- 
light radius of their galaxies was 0.3", so star/galaxy sepa- 
ration based upon morphology alone would be very difficult. 
They ascribe the excess in their number counts to inaccurate 
star/galaxy separation in the ground-based data. However, 
as noted above, our star/galaxy separation is based upon 
morphology and colour, and we do not find a large popu- 
lation of compact objects with the colours of galaxies. In 
addition, the CFRS collaboration obtained spectroscopy of 
one sixth of the objects in their sample, without regard to 
morphology, and also did not find a large population of com- 
pact galaxies. The excess in the MDS counts is most likely 
due to errors in the photometry of the WF/PC images in- 
troduced in the deconvolution process. 

The V— band galaxy counts are presented in Table [l] and 
Figure ^. We used the Kitt Peak standard V filter, corrected 
to the Johnson V standard stars of Landolt, (1983; 1992). 
There are only two other measurements of the V— band num- 
ber counts, by Casertano et al. (1995) (MDS), and by Driver 
et al. (1994). V555 is approximately equal to Johnson V for 
the galaxies in these samples, so we have applied no correc- 
tion. The MDS V counts again show an excess over the faint 
end of the counts presented here. 

The B— band galaxy counts are presented in Table B 
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Figure 3. The V— band number counts. The models plotted are 
as in Figure |lj. 

and Figure ^. We used the Kitt Peak standard B filter, cor- 
rected to the Johnson B standard stars of Landolt, (1983; 
1992). The B counts have been measured in many surveys, 
but this is the first time that B counts at B < 20 have been 
obtained with a CCD camera, rather than with non-linear 
photographic plates. While our area is much smaller than 
other surveys, and our statistical error is higher, our counts 
are far less likely to suffer from systematic effects in the pho- 
tometry. In Figure ^ we have plotted bj counts, converted 
to the Johnson B— band, from the APM survey of Maddox 
et al. (1990) (hereafter APM), and from the MAMA survey 
of Bertin & Dennefeld (1996). 

The APM counts show the steep slope at the bright 
end mentioned earlier. These authors interpretated the data 
as revealing a large (and unexpected) amount of luminos- 
ity evolution at low redshifts (2 < 0.1). Other workers have 
attributed this steep slope to a local underdensity of galax- 
ies (Shanks 1990; Metcalfe et al. 1991), to a selection effect 
against low surface brightness galaxies (McGaugh 1994; Fer- 
guson & McGaugh 1995), or to systematics in the photom- 
etry (Metcalfe et al 1995). The slope measured in the APM 
data at 16 < B < 19 is 0.59, while a linear fit to our data in 
this same range of magnitudes gives a slope of 0.50 ± 0.03. 
Our measured slope is consistent with passive evolution 
models, which have slopes of 0.52 and 0.51 for go = 0.5 
and respectively, and agrees with that of Bertin & Den- 
nefeld (1996), who surveyed 145 square degrees using indi- 
vidually calibrated Schmidt plates. We have normalized the 
flat-universe passive-evolution model with a least-squares 
fit to our data. This normalization is equivalent to using 
a Schechter luminosity function with b*j — —19.50 + 5log(h) 
and a — —0.97, as measured by Loveday et al. (1992), but 
with <))* — 2.02 x 10 -2 h 3 Mpc~ 3 , a normalization that is 
a factor of 1.44 times higher than that measured in the 
Stromlo-APM survey (Loveday et al. 1992.) High normal- 
ization models have been proposed to reduce the excess of 
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Figure 4. The B— band number counts. The models plotted are 
as in Figure All counts have been converted to Johnson B as 
discussed in the text. Our counts show a shallower slope than 
the APM counts of Maddox et al. (1990), and are consistent with 
passive evolution models. 



faint blue galaxies that has been seen in deep photometric 
surveys (for a review, see Metcalfe et al. 1996), and to fit 
the WFPC2 Medium Deep Survey results (Glazebrook et al. 
1995; Driver et al. 1995). The sensitivity of our survey to low 
surface brightness galaxies, and the effects of clustering on 
the errors of the number counts will be discussed elsewhere 
(Gardner et al., in preparation). 



4 CONCLUSIONS 

We have presented bright galaxy number counts measured 
with linear detectors in the B, V, I, and K bands in two 
fields totalling nearly 10 square degrees. All of our mea- 
surements are consistent with passive evolution models. Our 
counts do not exhibit the steep slope measured in other sur- 
veys, either in the K— or B— bands, and so do not support 
earlier interpretations that required a large amount of lumi- 
nosity evolution at low redshift. We also do not find evidence 
of a large underdensity in the local universe, unless it is a 
phenomenon occuring exclusively in the South Galactic Pole 
region. (Both of our fields are North of the Galactic Plane.) 
Our data are consistent with the conclusions of Metcalfe et al 
(1995) that the steep slope measured previously in the bright 
B— band number counts is most likely due to systematic er- 
rors in the non- linear photometry of photographic plates. 
Our B— band counts support the high normalization mod- 
els, based upon a local <f>* approximately 1.4 times higher 
than that measured by Loveday et al. (1992). Our if— band 
counts are consistent with previous measurements, and have 
reached a fair sample of the universe in the region where evo- 
lution and K-corrections are well understood. 
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